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INTRODUCTION 

Over  the  past  two  decades,  techniques  for  computing  complex  flows  have 
become  increasingly  more  sophisticated.  Steger  [1978],  Thompson  [1980],  Rubin 
and  Khosla  [1981,  1982],  Beam  and  Warming  [1982],  and  Briley  and  McDonald 
[1984]  have  computed  viscous  flows  at  reasonable  Reynolds  numbers;  Davis  and 
Werle  [1981]  and  Johnston  and  Sockol  [1984]  have  studied  viscid-inviscid 
interaction;  and  Edwards  and  Carter  [1985]  and  Melnik  and  Brook  [1985]  have 
computed  through  separated  regions.  Further,  all  computations  may  now  involve 
complex  turbulence  models,  such  as  the  models  by  Bradshaw,  Ferriss,  and  Atwell 
[1967]  and  Launder,  Reece,  and  Rodi  [1975].  It  is  desirable  that  these 
techniques  find  their  way  into  the  design  process.  The  design  process  of 
particular  interest  here  is  the  turbomachinery  design  process.  These 
numerical  techniques  are  capable  of  very  detailed  predictions  (for  example, 
boundary  layer  profiles),  but  to  be  used  with  confidence,  they  should  be 
tested  against  very  detailed  experimental  data  under  typical  flow  conditions. 
As  turbomachinery  testing  has  generally  been  concerned  with  overall 
turbomachinery  performance  rather  than  with  the  details  of  the  flow  field, 
such  data  are  lacking. 

In  order  to  provide  some  of  the  needed  data,  we  used  a  one-component 
laser  Doppler  velocimeter  (LDV)  to  measure  the  two-dimensional  periodic  flow 
field  about  a  double  circular  arc,  compressor  blade  in  cascade.  Eleven 
boundary  layer  profiles  were  taken  on  both  the  pressure  and  suction  surfaces 
of  the  blade;  two  profiles  were  taken  in  the  near  wake.  All  measurements  were 
made  at  a  chord  Reynolds  number  (Rec)  of  500,000  (±1%)  and  an  incidence  angle 
of  5  degrees  (that  is,  the  stagnation  point  is  on  the  pressure  surface).  The 
turbulence  intensity  in  the  incident  flow  was  0.18%.  With  an  incidence  angle 
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of  5  degrees,  the  pressure  surface  exhibits  a  large  region  of  laminar  flow  (up 
to  roughly  60%  chord);  transition  on  the  pressure  surface  appears  to  be 
incomplete.  The  suction  surface  profiles  appear  to  separate  both  at  the 
leading  edge  and  again  somewhat  beyond  midchord;  the  leading  edge  separation 
apparently  reattaches  by  2.6%  chord.  Using  the  terminology  of  Simpson,  Chew, 
and  Shivaprasad  (1981],  we  found  incipient  detachment  to  occur  at  60%  chord  on 
the  suction  surface  and  transitory  detachment  to  occur  at  83%  chord.  Inlet 
and  outlet  five-hole  probe  measurements  and  blade  static-pressure  measurements 
supplement  the  blade  boundary  layer  profiles.  Surface  flow  visualization, 
through  sublimation,  compliments  the  transition  and  separation  region  data. 
Before  describing  the  experiment  and  its  results,  we  will  describe  previous 
detailed  measurements  both  in  rotating  systems  and  in  cascades. 

II.  HISTORICAL  BACKGROUND 

Several  researchers  have  attempted  to  measure  boundary  layers  on 
turbomachine  blades.  Evans  (1978]  measured  boundary  layers  at  four  chord  (c) 
locations  at  midspan  on  the  suction  surface  of  a  stator  blade  (c  =  305.0  mm 
and  Rec  =  500,000).  The  hot-wire  measurements  were  made  in  an  axial-flow 
compressor  at  three  time-mean  incidence  angles  on  a  row  of  stationary  blades 
preceded  by  a  row  of  rotating  blades.  The  stator  blades  cut  the  wakes  of  the 
rotor  blades  and  the  rotor  wake  segments  are  subsequently  transported  through 
the  stator  passages.  Since  the  wake  segments  involve  low  velocity  fluid,  the 
boundary  layer  is  subject  to  a  periodically  varying  freestream,  and  the  blade 
is  subject  to  a  periodically  varying  incidence  angle.  Ensemble-averaged 
velocity  profiles  eliminate  the  random  unsteadiness  caused  by  turbulence. 
However,  the  periodic  unsteadiness  is  preserved.  The  ensemble-averaged 
velocity  profiles  at  30%  and  50%  chord  show  that  the  boundary  layers  alternate 
between  laminar  and  turbulent  because  of  the  unsteady  flow.  As  a  result,  the 


time-mean  velocity  profiles  exhibited  a  larger  boundary  layer  growth  than  was 
expected. 

Anand  and  Lakshminarayana  [1978]  measured  boundary  layers  on  the  rotor 
blades  of  a  rocket  pump  inducer  using  a  three-sensor  hot-wire  probe  rotating 
with  the  blades.  Because  of  imbalances  in  the  radial  pressure  force  and 
inertia  forces  in  the  blade  boundary  layer,  outward  radial  velocities  develop 
in  a  rotor  blade  boundary  layer  and  inward  radial  velocities  develop  in  a 
stator  blade  boundary  layer.  Anand  and  Lakshminarayana  [1978]  measured  a 
significant  outward  radial  component  in  the  boundary  layer  velocity  and  found 
that  this  radial  migration  strongly  influenced  the  chordwise  velocity 
profiles . 

Other  experimenters  have  measured  the  boundary  layers  on  rotor  blades  of 
axial-flow  fans.  Toyokura,  Kurokawa,  and  Kimoto  [1982]  used  rotating 
three-hole  cobra  probes  to  measure  the  three-dimensional  boundary  layers  at 
six  radial  sections  (c  =  80.0  mm  to  199.5  mm  and  Rec  *  300,000  to  500,000). 

The  outward  radial  flow  seemed  to  retard  the  predicted  regions  of  transition 
and  separation.  Lakshminarayana,  Govindan,  and  Hah  [1982]  used  rotating 
miniature  "x" -configuration  hot-wire  probes  for  boundary  layer  measurements  at 
five  radial  locations  (c  *  152.4  mm  and  Rec  =  280,000  at  midchord).  Rotating 
miniature  "x” -configuration  hot-wire  probes  were  also  used  by  Pouagare, 

Galmes,  and  Lakshminarayana  [1985]  for  numerous  boundary  layer  measurements  (c 
*  123.9  mm  to  154.1  mm).  In  no  case  could  the  velocity  profiles  be  very  well 
resolved. 

Walker  [1982]  made  measurements  similar  to  those  of  Evans  [1978].  Walker 
[1982]  used  a  hot  wire  in  an  axial-flow  compressor  to  measure  boundary  layers 
on  a  stator  (c  =*  76.0  mm  and  Rec  =  30,000  to  200,000)  downstream  of  both  a  row 
of  rotating  blades  and  a  row  of  inlet  guide  vanes.  He  tried  to  correct  for 


wall  proximity  using  a  method  outlined  by  Wills  [1962],  but  still  had 
difficulties  matching  the  law  of  the  wall.  Low  Reynolds  number,  large  adverse 
streamwise  pressure  gradients,  and  rapidly  changing  boundary  conditions  (due 
to  the  periodic  unsteadiness)  were  given  as  reasons  for  the  absence  of  a 
logarithmic  region. 

Hodson  [1983]  made  hot-wire  measurements  in  an  axial-flow  turbine  at 
midspan  on  a  rotor  blade  (c  =  114.5  mm  and  Rec  =  315,000)  downstream  of  a 
stator  row.  Once  again,  the  periodic  unsteadiness  seemed  to  cause  the 
boundary  layer  characteristics  to  vary  between  laminar  and  turbulent  at  some 
chordwise  locations.  Profile  losses  were  larger  than  expected  and  this  too 
was  attributed  to  the  unsteady  flow. 

Boundary  layer  measurements  on  turbomachine  blades  have  yet  to  produce 
velocity  profiles  with  enough  detail  and  precision  to  compare  with  viscous 
computational  codes.  Therefore,  our  understanding  of  the  physical  nature  of 
these  complex,  unsteady,  three-dimensional  boundary  layers  is  far  from 
complete.  These  flows  are  characterized  by  high  turbulence  levels,  as  well  as 
by  periodic  unsteadiness  caused  by  the  interaction  between  stationary  and 
rotating  blade  rows.  "Blockage"  effects  exist  because  'of  the  development  of 
the  end-wall  boundary  layers  and  the  consequent  contraction  of  the  mainstream 
flow.  The  blade  boundary  layers  are  also  arfected  by  centrifugal  and  Coriolis 
forces  associated  with  both  the  swirl  and  the  blade  rotation.  The  complex 
blade  geometries  and  the  complex  flow  field,  including  secondary  flows,  tip 
leakage,  and  trailing  vorticity,  make  the  analysis  or  measurement  of 
turbomachine  blade  boundary  layers  very  difficult. 

Because  of  these  difficulties,  many  researchers  have  sought  a  simple 
geometry  which  retains  some  of  the  physics  of  the  flow  in  which  to  make  their 
boundary  layer  measurements.  A  model  that  has  proven  effective  in  other  areas 


5 


of  turbomachinery  is  the  periodic,  two-dimensional  row  of  airfoils,  commonly 
referred  to  as  a  cascade.  Properly  realized,  a  cascade  should  eliminate  all 
of  the  complexities  of  the  turbomachine  except  blade  curvature,  secondary 
flow,  and  the  effect  of  freestream  turbulence. 

The  first  attempt  to  measure  boundary  layers  on  cascade  blades  was  made 
by  Peterson  [1958].  He  used  a  three-hole  cobra  probe  to  measure  blade 
boundary  layers  in  a  compressor  cascade  (c  =  123.8  mm  and  Rec  =  300,000). 
Several  boundary  layers  were  measured  on  both  the  suction  and  pressure 
surfaces  for  three  different  incidence  angles.  In  an  attempt  to  better  model 
an  actual  turbomachine,  Peterson  [1958]  simulated  the  added  diffusion  caused 
by  the  radial  distribution  of  axial  velocity  and  the  consequent  streamline 
deviation.  This  added  diffusion  was  created  in  the  cascade  by  placing  a 
perforated  metal  screen  downstream  of  the  blades.  The  measurements  were  not 
taken  in  the  freestream,  and  therefore,  the  normal  pressure  gradient  remains 
unknown.  This  lack  of  information  leads  to  a  problem  in  computing  the  edge 
velocity  (which  was  probably  inferred  from  the  blade  static-pressure 
distribution).  The  data  are  quite  scattered,  especially  in  the  regions  near 
separation,  where  the  cobra  probe  fails.  Peterson  [1958]  reported  no 
significant  differences  between  measurements  with  and  without  the  added 
diffusion . 

Pollard  and  Gostelow  [1967]  measured  boundary  layers  on  compressor 
cascade  blades  (c  =  152.4  mm  and  Rec  =  200,000)  with  a  Pitot  tube  to  examine 
the  effect  of  leading  edge  roughness.  Three  suction  surface  boundary  layers 
and  one  pressure  surface  boundary  layer  were  measured  on  both  blades  with  a 
smooth  leading  edge  and  with  a  leading  edge  roughened  with  polythene  spheres. 
With  no  roughness,  laminar  separation  occurred  before  transition.  It  appears 
that  separation  of  the  turbulent  boundary  layers  near  the  trailing  edge  would 
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|  be  more  likely  to  occur  when  leading  edge  roughness  is  present.  These  results 

j  agree  with  the  Preston  tube  skin-friction  measurements  of  Pollard  and  Gostelow 

i 

|  [1967].  The  detail  and  precision  of  the  data  are  unknown  since  no  actual  data 

|  points  are  given. 

|  Evans  [1981]  used  compressor  cascade  blades  (c  ■  304.8  mm  and 

' 

Rec  *  500,000)  for  boundary  layer  measurements  taken  with  a  hot-wire 
anemometer  probe.  A  problem  with  this  data  is  that  the  blade  boundary  layers 
were  tripped  with  a  wire  at  10%  chord.  Evans  [1971]  argued  that  a  turbulent 
boundary  layer  over  most  of  the  blade  would  better  represent  the  high 
turbulence  and  unsteadiness  levels  usually  encountered  in  a  turbomachine. 
Instead,  the  artificially  induced  boundary  layer  development  is  rather 
misleading. 

Problems  of  contamination,  corrosion,  erosion,  and  deposition  have  led  to 
two  investigations  that  dealt  with  the  effects  of  surface  roughness  on  cascade 
blade  boundary  layers.  Bammert  and  Milsch  [1972]  measured  blade  boundary 
layers  in  a  compressor  cascade  (c  =  180.0  cm  and  Rec  =  430,000)  with  four 
different  blade  profiles  to  parametrically  change  the  pitch,  camber,  and 
thickness  of  the  blades.  They  used  emery  powder  to  develop  the  five  roughness 
grades  to  be  tested.  A  turbine  cascade  (c  =  17.50  cm  and  Rec  =  560,000)  was 
used  for  the  blade  boundary  layer  measurements  of  Bammert  and  Sandstede  [1980] 
where  four  roughness  grades  were  tested.  Both  investigations  used  a  flattened 
Pitot  tube  which  allowed  measurements  to  be  taken  very  close  to  the  blade 
surface.  Bammert  and  Sandstede  [1980]  used  a  hot-wire  anemometer  to  confirm 


the  measurements  made  with  the  flattened  Pitot  tube.  The  studies  showed  that 
increasing  surface  roughness  led  to  increases  in  both  momentum  thickness  and 
skin  friction  and  a  forward  shift  of  the  regions  of  transition  and  separation 


Meauze  [1979]  used  a  transonic  compressor  cascade  (c  =  9.49  cm, 

Rec  ■  1,660,000  for  Mj  *  0.70,  and  Rec  *  2,120,000  for  Mj  ■  0.85)  for  blade 
boundary  layer  measurements.  He  used  total-pressure  probes  to  measure  suction 
surface  boundary  layers  for  two  inlet  Mach  numbers  and  four  incidence  angles. 
The  flow  recompression  on  the  highly  cambered  blades  results  in  laminar  flow 
separation  and  subsequent  turbulent  reattachment.  The  thin  laminar  boundary 
layers  upstream  of  this  separation  "bubble"  were  difficult  to  measure,  so  that 
the  total-pressure  profiles  were  only  measured  for  turbulent  boundary  layers. 

The  trailing  edge  boundary  layers  on  both  the  suction  and  pressure 
surfaces  of  a  compressor  cascade  blade  (c  =  203.2  mm  and  Rec  =  478,000)  were 
measured  with  a  hot-film  probe  by  Hobbs,  Wagner,  Dannenhof f er ,  and  Dring 
[1980].  The  two  profiles  are  very  detailed  and  precise;  they  show  a  nearly 
separated  profile  on  the  suction  surface,  and  also  a  pressure  surface  profile 
typical  of  favorable  pressure  gradients. 

Recently,  Hodson  [1983]  has  presented  blade  boundary  layer  data  measured 
in  a  turbine  cascade  (c  =  114.5  mm  and  Rec  =  315,000)  with  a  hot-wire  probe. 
Although  no  data  points  are  reported,  the  velocity  profiles  on  the  suction 
surface  show  laminar  flow  until  78%  chord  followed  by  laminar  separation  and 
no  reattachment. 

All  prior  investigations  have  used  hot-wire  or  pressure  probes  to  make 
boundary  layer  measurements.  There  are  three  potential  problems:  first,  a 
periodic  twc-dimensional  cascade  flow  is  difficult  to  establish  and  the  probe 
may  distort  it;  second,  some  boundary  layers  are  likely  to  be  small  compared 
to  the  prole  dimensions  (some  may  be  small  relative  to  the  LDV  measurement 
volume,  see  Section  IV);  and  third,  separated  regions,  if  present,  cannot  be 
conveniently  studied.  For  these  reasons,  the  current  study  used  the  LDV 
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technique.  Before  discussing  the  LDV  results,  however,  we  will  describe  the 
cascade  facility  and  flow. 

III.  CASCADE  TESTS 

The  cascade  wind  tunnel  has  a  0.37  m  by  0.64  m  test  section  with  a 
maximum  air  speed  of  approximately  35  m/sec.  Tunnel  turbulence  control  is 
through  a  honeycomb  with  a  3.18  mm  cell  size,  several  settling  screens  in  the 
diffuser  section,  and  a  nine-to-one  contraction.  Over  the  speed  range  of 
24-35  m/sec,  the  tunnel  operated  with  a  freestream  turbulence  level  of  0.18% 
±10%  as  measured  by  a  hot-wire  anemomenter.  Figure  1  shows  a  schematic  of  the 
open  return  facility. 

The  blade  section  used  in  the  tests  is  a  compressor  blade  designed  at  the 
NASA  Lewis  Research  Center  (see  Sanger  [1980]).  The  blade  section  is  a  double 
circular  arc  blade  with  65  degrees  of  camber*  and  a  228.6  mm  chord  length. 

Both  the  leading  and  trailing  edges  were  machined  to  a  9.14  pm  radius.  The 
blade  aspect  ratio  is  1.61.  The  five  cascade  blades  were  made  of  aluminum  and 
were  anodized  black  to  minimize  laser  reflections.  Camber  line  and  thickness 
relationships  necessary  to  construct  the  blades  (or  for  computation)  are  given 
in  the  Appendix.  To  insure  proper  alignment,  the  blades  were  carefully 
positioned  in  two  inserts  (aluminum  and  plexiglass**  (on  the  optics  side)), 
which  were  in  turn  mounted  to  the  plywood  walls  of  the  cascade  test  section. 
The  test  section  is  shown  in  Figure  2.  For  the  case  to  be  reported,  the 
cascade  had  a  solidity  of  2.14.  The  stagger  angle  was  20.5  degrees.  The 
important  cascade  and  flow  angles  are  defined  in  Figure  2. 


A  design  condition  was  that  the  blade  show  some  trailing  edge  separation  of 
the  suction  surface  turbulent  boundary  layer  at  zero  incidence  angle.  With 
the  available  cascade  geon*.  r  y,  this  leads  to  a  large  camber  angle. 

*A  glass  insert  was  later  placed  within  the  plexiglass  to  improve  the  I.DV 
signal. 
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As  current  computer  codes  assume  a  two-dimensional  periodic  cascade  flow. 


one  must  take  data  in  such  a  flow  field  for  it  to  be  of  use.  Two- 


dimensionality,  of  course,  implies  that  the  velocities  and  angles  of  the  flow 


are  substantially  the  same  in  spanwise  planes,  while  periodicity  supposes  that 


velocities  and  flow  angles  show  only  minimal  variations  from  blade  passage  to 


blade  passage,  both  upstream  and  downstream  of  the  blade  row.  For  the 


five-bladed  cascade  used  here,  periodicity  was  taken  to  mean  periodicity  over 


three  blade  passages  centered  at  the  minimum  velocity  point  of  the  middle 


blade  wake.  In  order  to  satisfy  the  condition  of  continuity  for  a 


two-dimensional,  incompressible  cascade  flow,  the  axial  velocity  must  be  held 


constant  throughout.  From  a  practical  standpoint  then,  flow  periodicity  and 


two-dimensionality  could  be  determined,  for  a  uniform  inlet  flow,  by  examining 


the  axial-velocity  ratio  and  flow  angles  determined  at  the  exit  plane.  This 


was  quite  useful  as  these  outlet  measurements  could  be  made  simply  and 


quickly,  thus  allowing  us  to  check  the  cascade  flow  daily. 


Two-dimensionality  and  periodicity  are  normally  controlled  in  cascades  by 


using  continuous  suction  over  the  entire  blade  pack.*  Continuous  suction  was 


not  possible  in  the  current  experiment  because  of  the  need  for  laser  access. 


Alternate  flow  control  was  examined  in  some  detail.  Returning  to  Figure  2,  we 


note  that  there  are  many  potentially  useful  flow  controls.  That  is,  it  is 


possible  in  principle,  to  control  the  flow  by  adjusting  the  position  of  the 


lower  false  blade,  the  upper  false  blade,  the  variable  diffuser,  and  the 


tailboards  (as  well  as  the  relative  position  of  the  tailboards),  or  by 


adjusting  the  magnitude  (and  distribution)  of  the  top  suction,  side  suction, 


and  lower  and  upper  channel  suction.  In  practice,  the  lower  and  upper  false 


Note  that  the  control  here  seems  to  be  over  the  growth  of  channel  corner 
disturbances.  An  attempt  to  use  a  slightly  divergent  wall,  commonly  used  to 
compensate  for  boundary  layer  growth,  failed. 
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blades  were  each  set  at  nominally  one  blade  spacing  from  blades  1  and  5, 
respectively.  The  diffuser  was  set  to  minimize  the  flow  angle  at  the  splitter 
plate.  Lower  channel  suction  was  not  required.  Top  and  upper  channel  suction 
were  provided  through  the  same  5  hp  blower,  which  was  run  at  full  power.  A 
baffle  system  was  used  to  adjust  the  relative  amounts  of  suction  provided  at 
the  top  and  upper  channel,  and  the  baffling  along  with  slight  adjustments  to 
the  upper  false  blade  position  were  used  to  insure  a  horizontal  flow  at  the 
upper  false  blade  leading  edge.  The  tailboards  were  most  useful  in 
controlling  the  relative  exit  angles  of  the  flow;  that  is,  they  could  be 
adjusted  so  that  the  exit  angle  across  the  cascade  in  a  blade-to-blade 
direction  was  constant  outside  of  the  wake  regions.  The  periodicity,  however, 
was  found  to  be  most  influenced  by  the  amount  and  distribution  of  the  sidewall 
suction.  Side  suction  was  provided  by  a  10  hp  centrifugal  blower  operated  at 
full  power.  Side  suction  distribution  was  controlled  by  a  complex  baffling 
system.  Six  suction  ducts  were  located  at  a  half  blade  spacing  on  either  side 
of  each  of  the  blades.  Each  of  the  individual  ducts  had  a  separate  baffle 
control,  and  was  adjusted  by  the  simple  but  tedious  procedure  of  changing  a 
baffle  position  and  then  examining  the  resulting  outlet  flow.  Presumably,  the 
control  of  side  suction  distribution  controlled  the  size  of  the  sidewall 
boundary  layer  at  its  intersection  with  the  blade  pack  leading  edge  line — in  a 
sense  controlling  the  virtual  origin  of  the  corner  disturbances.  Control  of 
sidewall  suction  distribution  then  implied  a  control  of  the  "blockages,” 
caused  by  corner  disturbance  contamination,  of  each  blade  passage  individually 
and  hence  control  of  the  individual  blade  angles  of  attack.  Once  set,  the 
stability  of  the  periodicity  obtained  on  a  day-to-day  basis  was  excellent. 

The  two-dimensionality  and  periodicity  of  the  cascade  flow  were 
determined  by  conducting  measurements  of  the  cascade  inlet  and  outlet  velocity 


profiles  using  five-hole  probes.  These  probes  are  capable  of  resolving  the 
three  components  of  velocity  as  well  as  the  relative  yaw  and  pitch  angles  of 
the  flow.  Treaster  and  Yocum  [1979]  give  a  complete  description  of  the 
five-hole  probes  employed  in  the  study.  The  probes  were  calibrated  at  a  speed 
of  30.5  m/sec  in  an  open  jet  air  facility  over  the  range  of  pitch  and  yaw 
angles  of  +30  degrees  to  -30  degrees.  Reynolds  number  effects  on  the  probes 
have  been  shown  to  be  small  (see  Treaster  and  Yocum  [1979]). 

All  inlet  and  outlet  five-hole  probe  surveys  were  referenced  to  a 
Pitot-static  probe  located  25.4  mm  upstream  of  the  blade  pack  leading  edge 
line.  The  probe  protruded  approximately  150  mm  into  the  flow — well  outside  of 
the  sidewall  boundary  layer.  The  probe  was  located  along  a  parallel  to  the 
leading  edge  line  at  a  position  between  blades  3  and  4  at  which  the  inlet 
velocity  was  roughly  equal  to  the  average  of  the  measured  inlet  velocities. 

Cascade  inlet  flow  profiles  were  documented  by  five-hole  probe 
measurements  approximately  38  mm  upstream  of  and  parallel  to  the  leading  edge 
line.  Outlet  flow  profiles  were  measured  parallel  to  the  trailing  edge  line 
and  one-half  chord  downstream  of  this  line.  A  nearly  real-time  data 
acquisition/reduction  system  was  used  for  the  velocity  measurements.  In  each 
case  the  five-hole  probe  and  Pitot-static  probe  data  were  sent  to  seven 
separate  pressure  transducers.  These  seven  pressure  signals  as  well  as  a  test 
section  temperature  signal  were  scanned  by  a  raultiplexer/scanner ,  smoothed  on 
a  multimeter  through  a  100  cycle  (1  2/3  sec)  integration,  and  sent  for 
reduction  to  a  VAX  11/782.  Velocities,  velocity  ratios,  and  flow  angles  could 
be  displayed  on  a  video  terminal,  written  on  a  line  printer,  or  plotted  on  a 
flat-bed  plotter. 

Inlet  velocity  surveys  were  made  after  the  outlet  flow  had  been 


determined  to  be  satisfactory.  Good  periodicity  and  spanwise  consistency  were 
apparent,  as  was  some  streamline  bending  induced  by  the  presence  of  the 
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blades.  An  average  Inlet  velocity  was  found  to  be  33.11  m/sec,  with  an 
average  incidence  angle  of  5  degrees.  Chord  Reynolds  number,  based  on  the 
inlet  velocity,  was  500,000  with  an  observed  ±1%  variation  on  a  day-to-day 
basis.  Figure  3  shows  a  typical  outlet  flow  profile  and  the  equivalent 
turning  angles.  The  periodicity  of  the  flow  is  clearly  excellent.  The 
axial-velocity  ratio  is  determined  by  averaging  the  local  axial-velocity  ratio 
over  three  blade  passages,  centered  at  the  minimum  velocity-ratio  point  of  the 
center  blade  wake.  In  the  calculation,  the  inlet  flow  is  assumed  to  be 
spatially  constant  at  its  average  value  (the  Pitot-static  tube  reading).  The 
average  axial-velocity  ratio  was  found  to  be  1.0;  on  a  day-to-day  basis  the 
variation  was  within  ±3%.  The  flow  turning  angle  averaged  across  the  three 
center  blade  passages  was  54  degrees.  Figure  4  shows  all  of  the  blade 
geometry  and  the  inlet  and  outlet  flow  measurements. 

Measurements  have  been  taken  to  help  quantify  the  losses  in  total 
pressure  across  the  cascade.  Additional  quantities  can  be  computed  to  compare 
with  design  limits  on  diffusion  rate  and  static-pressure  rise  within  the 
cascade.  The  difference  between  the  blade-passage-averaged  flow  angle  in  the 
outlet  flow  and  the  exit  blade  "metal"  angle  is  the  deviation  angle.  For  the 
current  study,  the  deviation  angle  was  measured  to  be  16  degrees,  which  is 
very  large.*  Non-dimensionalizing  by  the  inlet  dynamic  pressure,  the 
blade-passage-averaged  total-pressure  loss  coefficient  was  0.151  and  the 
blade-passage-averaged  static-pressure  rise  coefficient  was  0.463.  Note  that 
the  static-pressure  rise  coefficient  was  probably  affected  by  the  positioning 
of  the  tailboards.  An  equation  for  the  total-pressure  loss  coefficient  was 
developed  by  Lieblein  and  Roudebush  [19561  where 


*A1 though  large,  a  16  degree  deviation  angle  is  not  unexpected,  as  the  design 
condition  called  for  a  zero  incidence  angle.  Experiments  at  lower  incidence 
angles  are  planned. 
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A  value  of  0.172  for  u  can  be  calculated  from  this  equation  using  the  flow 
parameters  measured  in  the  outlet  flow.  The  loss  in  total  pressure  across  the 
cascade  is  related  to  the  amount  the  flow  is  diffused  through  the  blade 
passage.  Lieblein,  Schwenk,  and  Broderick  [1953]  derived  a  diffusion  factor 
where 

V2  “  ^02  cos  sin  6l  -  cos  tan  62 

^  *  V1  +  2 oVj  *  cos  62  +  2o 

Designers  normally  place  a  limit  of  0.6  on  D  when  designing  a  blade  row.  This 
limit  is  based  on  a  large  number  of  cascade  performance  tests.  Values  of  D 
greater  than  0.6  result  in  large  increases  in  total-pressure  loss  because  the 
larger  amount  of  diffusion  causes  the  blade  boundary  layers  to  separate. 

Using  a  blade-passage-averaged  value  of  the  outlet  flow  angle  allows  a  value 
of  0.658  for  D  to  be  computed  for  the  current  study.  This  value  of  D 
indicates  a  risk  of  separation. 

Ideally,  blade  static  pressure  would  be  measured  on  the  center  blade  of 
the  cascade — the  one  intended  for  LDV  measurements — but  the  two  types  of 
measurements  had  somewhat  conflicting  requirements.  That  is,  an 
aerodynamical ly  smooth  (0.8  pm  estimated  surface  roughness)  surface  was 
desired  for  the  LDV  surveys,  while  the  conveying  tubes  required  for  the 
pressure  measurements  inevitably  led  to  a  somewhat  roughened  surface.  To  work 
around  this  problem,  we  instrumented  the  suction  surface  of  the  upper  blade,  4 
(see  Figure  2),  and  the  pressure  surface  of  the  lower  blade,  2,  with  24 
pressure  taps;  we  instrumented  the  center  blade,  3,  pressure  and  suction 
surfaces  with  6  and  7  taps,  respectively.  Since  the  flow  was  periodic,  the 
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pressure  distribution  could  be  obtained  from  blades  2  and  4  as  well  as  from 
blade  3.  After  checking  this  by  comparing  the  results  of  the  pressure 
distributions  from  blades  2  and  4  against  the  data  from  blade  3,  (the 
agreement  was  excellent),  we  interchanged  blade  3  with  the  uninstrumented 
blade,  5,  for  the  LDV  surveys.  Data  acquisition  and  reduction,  for  the 
static-pressure  distribution,  were  similiar  to  that  described  for  the 
five-hole  probe  data  with  the  exception  that  a  scanner  valve  was  used  to 
switch  the  pressure  data,  hole-by-hole,  to  a  single  transducer  during  data 
acquisition. 

The  blade  static-pressure  distribution  appears  in  Figure  5.  Integrating 
this  distribution  gives  a  lift  coefficient  of  0.952.  The  5  degree  incidence 
angle  dramatically  affects  the  pressure  distribution.  The  large  favorable 
gradient  on  the  pressure  surface  suggests  that  the  boundary  layers  near  the 
leading  edge  will  be  laminar;  transition  should  be  looked  for  on  the  pressure 
surface.  The  unfavorable  gradient  at  the  leading  edge  of  the  suction  surface 
implies  a  leading  edge  separation.  The  rather  flat  pressure  profile  near  the 
trailing  edge  of  the  suction  surface  suggests  that  the  flow  may  be  separated 
there.  The  rapid  and  continuous  changes  of  pressure  on  both  the  pressure  and 
suction  surfaces  offer  little  hope  of  finding  equilibrium  boundary  layers.  In 
addition,  the  inviscid  velocity  field  within  the  blade  passage  will  be  under 
the  influence  of  a  normal  pressure  gradient,  and  one  cannot  anticipate  a 
constant  freestream  velocity  region. 

Following  the  sublimation  method  used  by  Holmes  and  Obara  [1982],  we 
used  surface  flow  visualization  to  help  determine  the  location  of  the 
anticipated  trailing  edge  separation  on  the  suction  surface,  and  the 
transition  on  the  pressure  surface.  An  air  brush  was  used  to  coat  the  blade 
with  a  mixture  of  napthalene  and  acetone  in  a  1:8  solid-to-solvent  volume 
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ratio.  Preliminary  tests  determined  that  a  run  time  of  nine  minutes  was 
sufficient  to  set  up  the  visualization  pattern.  For  each  run  the  center  blade 
was  removed,  the  napthalene/acetone  mixture  was  applied,  and  large  particles 
were  dusted  off.  The  blade  was  then  replaced  in  the  cascade  and  removed  after 
the  nine  minute  run  time,  allowing  the  sublimation  pattern  to  be  photographed. 
A  sufficient  number  of  tests  were  taken  so  that  a  meaningful  mean  and 
Student's  t  test  deviation  could  be  obtained.  Evidence  of  two-dimensionality 
was  much  more  apparent  in  the  suction  surface  separation  pattern  then  in  the 
pressure  surface  transition  pattern.  This  is  perhaps  an  indication  of  the 
importance  of  the  local  surface  roughness  in  determining  the  transition  point 
for  the  very  thin  boundary  layers  encountered.  With  95%  confidence, 
separation  was  found  to  occur  at  65.6%  chord  on  the  suction  surface  with  a 
deviation  of  ±3.5%  chord.  Transition  was  found  on  the  pressure  surface  to  be 
at  64.2%  ±  3.9%  chord  to  the  same  level  of  confidence. 

IV.  THE  LASER  DOPPLER  VELOCIMETER 

All  blade  boundary  layer  and  near-wake  measurements  were  made  using  a 
single-component  LDV.  For  all  the  LDV  measurements,  a  specially  designed 
traversing  mechanism  was  used  which  matched  the  arc  of  motion  of  an  optics 
cradle  to  that  of  the  blade  curvature  (two  arcs  were  employed,  one  for  each  of 
the  pressure  and  suction  surfaces).  All  measurements  were  made  in  the  plane 
of  the  local  blade  normal.  Translation  of  the  optics  cradle  normal  to  the 
blade  could  be  accomplished  in  step  intervals  as  small  as  0.0254  mm.  Prior  to 
the  LDV  measurements,  a  reference  distance  was  established  by  focusing  the  LDV 
control  volume  on  an  insert  which  securely  fit  over  the  center  measuring 
blade.  Narrow  lines  had  been  etched  on  the  insert  (on  arcs  matching  the  blade 
curvature)  to  be  known  distances  from  the  blade  surface.  Repeatability  in 
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establishing  a  measurement  reference  was  estimated  to  be  ±0.05  mm,  and  this 
uncertainty  is  probably  the  major  source  of  scatter  in  the  velocity  data. 

A  two-watt,  Spectra-Physics,  Argon-Ion  laser  was  used  for  the 
measurements.  Power  on  the  blue  line  (488  nm  wavelength)  ranged  between  0.5 
watts  and  0.7  watts  on  a  day-to-day  basis.  Standard  TSI  backscatter  optical 
components  were  used:  the  371.3  mm  focusing  lens  was  chosen  to  allow  the 
measurements  to  be  made  at  the  blade  midspan.  The  ellipsoidal  measurement 
volume  was  reduced  through  the  use  of  a  (2.71:1)  beam  expander;  the  predicted 
measurement  length  in  the  normal  to  the  blade  direction  was  37  pm.  While  this 
length  was  small  when  compared  to  the  length  scales  of  the  turbulent  boundary 
layers  on  the  suction  surface,  we  shall  show  that  it  is  roughly  half  the  size 
of  the  displacement  thickness  of  the  initial  laminar  profile  on  the  pressure 
surface.  Where  appropriate,  optical  shifting  at  5  MHz  was  employed.  Note 
that  to  measure  close  to  the  blade  surface  the  optical  cradle  was  tilted  at  an 
angle  of  roughly  1  degree.  Silicon  carbide  particles  having  a  mean  diameter 
of  1.5  pm  were  used  for  laser  seeding.  In  an  attempt  to  maintain  a  uniform 
distribution,  we  injected  the  silicon  carbide  particles  well  upstream  of  the 
measurement  station  at  the  flexible  coupling  (see  Figure  1).  The  particles 
were  suspended  in  a  "cloud  chamber,"  which  was  constructed  for  this  study,  and 
were  injected  into  the  tunnel  by  a  small  overpressure.  Gain,  on  the  counter 
processor,  was  kept  low,  and  particle  counts  averaged  only  20  or  so  particles 
per  second  (as  we  will  discuss,  however,  the  velocity  probability 
distributions  were  remarkably  clean). 

LDV  data  acquisition  and  reduction  was  accomplished  by  using  a  direct 
link  to  a  Vax  11/782  computer.  Software  allowed  selection  of  the  focusing 
lens  half  angle,  the  laser  wavelength,  the  frequency  shift,  the  minimum  number 
of  cycles  employed  in  the  calculation  (8  here),  and  the  number  of  particle 


counts  per  run.  Initial  output  was  in  the  form  of  a  velocity  histogram. 
Minimum  and  maximum  velocity  limits  could  be  set  by  two  cursors  to  eliminate 
obvious  noise  from  the  distribution.  Final  output  was  mean  velocity,  local 
turbulence  intensity,  and  the  percent  of  particle  counts  employed  in  the 
calculations.  For  some  of  the  profiles  measured  the  skewness  and  kurtosis  of 
the  distribution  were  also  calculated.  The  percentage  of  particle  counts 
employed  in  the  calculation  may  be  used  as  an  indicator  of  signal-to-noise 
ratio.  At  least  98%  of  the  total  particle  counts  were  used  for  measurement 
stations  in  the  boundary  layer;  at  least  95%  were  employed  for  points  in  the 
freestream  (the  difference  in  percentages  reflects  the  fact  that  fewer  overall 
points  were  used  at  the  freestream  locations). 

For  a  counter  processor,  employed  in  a  highly  turbulent  flow,  the 
calculation  of  mean  velocity  and  turbulence  intensity  may  not  be 
straightforward.  McLaughlin  and  Tiederman  [1973],  Hoessel  and  Rodi  [1977], 
Giel  and  Barnett  [1979],  Edwards  [1981],  Edwards  and  Jensen  [1983],  Johnson, 
Modarress,  and  Owen  [1984],  and  Stevenson,  Thompson,  and  Craig  [1984]  have  all 
discussed  the  question  of  velocity  bias  in  a  highly  turbulent  flow.  As 
pointed  out  first  by  McLaughlin  and  Tiederman  [1973],  the  problem  arises 
because  more  high-speed  particles  than  low-speed  particles  arrive  in  the 
measurement  volume  during  a  given  measurement  interval.  A  related  problem, 
termed  incomplete  signal  bias  by  Stevenson,  Thompson,  and  Craig  [1984],  can  be 
eliminated  by  employing  a  sufficiently  high  frequency  shift. 

McLaughlin  and  Tiederman  [1973]  describe  a  correction  for  the  phenomenon 
(with  uniform  seeding) — but  it  requires  complete  velocity  vector  Information. 

A  more  practical  one-dimensional  correction  is  also  given,  but  this  correction 
tends  to  overestimate  the  error  for  local  turbulent  intensities  >20%.  Edwards 
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[1980]  shows  that  the  biasing  error  can  be  made  negligibly  small  for  the  case 
of  a  saturated  data  handling  system  (taking  data  at  a  fixed  rate  set  by  the 
slowest  component  of  the  system)  by  taking  the  particle  density  equal  to  the 
seeding  rate  (this  assumes  a  validation  circuit  for  the  system).  Stevenson, 

Thompson,  and  Craig  [1984]  used  equal  time  interval  sampling  in  a  very  highly 
seeded  mixing  layer  as  a  bias-free  test  case.  As  pointed  out  by  Edwards  and 
Jenson  [1983],  however,  very  high  seeding  rates  may  open  the  door  to  other 
types  of  errors,  for  example,  by  reducing  the  actual  number  of  statistically 
independent  samples  used  to  form  the  velocity  statistics.  Moreover,  a  very 
high-seed  rate  may  be  very  difficult  to  achieve  in  precisely  those  regions  in 
which  the  bias  is  expected  to  be  high.  Often,  in  fact,  a  counter  processor  is 
chosen  over  a  tracker  processor  because  of  its  ability  to  act  at  very  low 
seeding  rates. 

The  bias  question  is  obviously  quite  complex,  and  a  consensus  opinion  on 
how  to  correct  data  is  still  lacking.  Some  issues,  such  as  non-uniform 
particle  seeding,  of  interest  particularly  in  air  flows,  have  not  yet  been  the 
subject  of  detailed  studies  (see  Hoessel  and  Rodi  [1977]  for  instance).  Giel 
and  Barnett  [1979]  conducted  an  experiment  favorable  to  obtaining  statistical 
bias,  but  no  consistent  bias  was  evident — thus  further  obscuring  the  bias 
question.  In  the  current  study,  we  employed  simple  arithmetic  averaging.  For 
many  of  the  boundary  layers  measured,  we  monitored  the  skewness  with  the  idea 
that  a  change  of  shape  from  the  classical  distribution  in  the  boundary  layer 
might  signal  significant  velocity  bias.  No  such  deviations  were  observed.  We  ^ 

note  that  both  Mclaughlin  and  Tiederman  [1973]  and  Johnson,  Modarress,  and 
Owen  [1984]  show  that  the  overestimate  of  the  mean  velocity  goes  roughly 
quadratically  with  turbulence  intensity,  being  5%  for  a  local  turbulence 
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intensity  near  20%,  and  being  12%  for  a  local  turbulence  intensity  near  35%.* 
These  numbers  should  be  borne  in  mind  not  only  when  evaluating  the  data 
presented  here,  but  also  when  evaluating  any  measurements  made  in  highly 
turbulent  flows.  The  mean  velocity  here  was  taken  as 

1  N 

u  I ,  Un 

n=l 


the  local  turbulence  intensity  was  taken  as 

1/2 


u'  1  fl  r  ,  v 

~  =  u  N  i,  (un  ~ 

L  n=1 


and  the  turbulence  intensity  was  taken  as  u'/Ue. 

Experience  has  shown  that  quite  satisfactory  repeatability  of  the  mean 
velocity  and  turbulence  intensity  can  be  guaranteed  in  boundary  layer  flows  by 
using  1000  particle  counts  in  regions  in  which  the  local  turbulence  intensity 
exceeds  5%.  In  regions  of  local  turbulence  intensity  of  less  than  5%  but  more 
than  2%,  500  points  are  used,  while  200  points  are  used  in  regions  of  less 
than  2%.  At  each  chord  position,  profiles  were  defined  by  statistically 
treating  the  data  for  six  individual  experiments.  Six  experiments  were  chosen 
as  the  statistics  found  from  six  experiments  showed  less  than  1%  scatter  in 
the  freestream  data.  Error  bands,  presented  on  the  LDV  data  plots,  represent 
95%  confidence  levels  as  determined  by  a  Student's  t  test. 

The  preliminary  data  analysis  is  automated  on  the  VAX  11/782  computer. 

The  effect  of  the  normal  pressure  gradient  on  the  boundary  layer  profiles  is 
accounted  for  first.  Details  of  the  technique  are  given  by  Zierke  and  Deutsch 
[1985]  and  basically  follow  the  approach  suggested  by  Mellor  and  Wood  [1971] 


*It  is  not  obvious  that  turbulence  intensity  is  the  only  relevant  variable. 
For  example,  the  values  of  the  higher  moments  of  the  velocity  probability 
distribution  are  no  doubt  also  of  importance. 


and  Ball,  Reid,  and  Schmidt  [1983].  Briefly  the  technique  assumes  that  the 
profile  may  be  represented  as 

u  "  umeas  “  uinv  + 

so  that  the  edge  velocity  (Ue)  can  be  determined  by  extrapolating  the  outer 
inviscid  flow  (u^nv)  to  the  wall  (where  u  =  umeas  =  0)  in  some  reasonable 
manner.  The  method  is  not  rigorous  in  its  definition  of  the  inviscid  region 
and  hence  in  the  manner  of  extrapolation.  Our  own  experience  with  the  22 
boundary  layers  measured  here,  however,  is  that  the  edge  velocity  is  quite 
insensitive  to  any  reasonable  choice  of  the  inviscid  region. 


V.  PRESSURE  SURFACE  BOUNDARY  LAYERS 


Boundary  layer  measurements  were  made  at  11  chord  locations  on  the 
pressure  surface  of  the  center  cascade  blade.  To  help  interpret  these 
velocity  profiles.  Figure  5  shows  the  measurement  locations  along  with  the 
pressure  distribution.  The  combination  of  continuously  changing  pressure  and 
moderate  surface  curvature  (0.002  <  |6/RC|  <  0.02)  signals  a  complicated 
non-equilibrium  flow  field.  At  the  leading  edge,  for  example,  the  large 
incidence  angle  (5  degrees)  results  in  a  strong  acceleration  which  promises  a 
region  of  laminar  flow.  In  the  region  from  8%  chord  to  62%  chord,  the  flow  is 
subjected  to  a  mildly  adverse  gradient  so  that  the  onset  of  transition  might 
be  expected  in  this  region.  The  subsequent  favorable  gradient,  however,  makes 
the  eventual  complete  transition  to  turbulence  problematic. 

The  measured  pressure  surface  boundary  layers  are  shown  in  Figure  6.*  The 
blade-to-blade  pressure  gradient  affects  the  inviscid  region  of  each  profile. 
This  pressure  gradient  varies  from  a  strong,  non-linear  gradient  near  the 
leading  edge,  where  the  streamlines  have  a  large  curvature,  to  a  nominally 


*Tabulated  data  will  be  supplied  upon  request. 


zero  gradient  near  the  trailing  edge.  As  previously  noted,  each  profile  was 
measured  six  times  and  the  symbols  in  Figure  5  represent  velocity  data 
averaged  over  the  six  tests.  The  error  bands  give  95%  confidence  levels  as 
determined  by  the  Student’s  t  test.  These  error  bands  are  quite  small, 
particularly  in  the  inviscid  regions.  The  exception  appears  at  2.7%  chord 
where  the  boundary  layer  is  so  small  that  the  LDV  could  only  nominally 
penetrate  the  layer,  and  correspondingly,  the  resolution  is  poor. 

The  boundary  layers  were  analyzed  using  methods  described  by  Zierke  and 
Deutsch  [1985].  The  influence  of  the  normal  pressure  gradient  was  first 
removed.  The  reconstructed  boundary  layer  data  were  then  compared  with  a 
Falkner-Skan  velocity  profile  (see  Falkner  and  Skan  [1931])  at  the  local 
streamwise  pressure  gradient.  For  the  velocity  profiles  measured  near  the 
trailing  edge,  an  attempt  to  fit  the  results  to  the  wall-wake  equation  of 
Coles  [1956]  was  also  made.  Finally,  integral  parameters  were  obtained  from 
both  a  smoothed  cubic  spline  fit  of  the  data  and  from  the  Falkner-Skan 
solutions.  The  velocity  profiles  are  replotted  non-dimensionally  in  Figure  7. 
In  spite  of  the  influence  of  both  curvature  and  changing  pressure  gradient  on 
the  flow  field,  the  Falkner-Skan  approximation  appears  to  reasonably  represent 
the  mean  velocity  profiles  through  about  57.2 %  chord.  At  68.0%  chord  and 
beyond,  there  is  an  increased  thickening  of  the  measured  profiles  relative  to 
the  Falkner-Skan  correlation  which  indicates  transitional  boundary  layers. 

Empirical  relationships  have  been  developed  for  the  prediction  of  the 
beginning  and  end  of  transition;  these  relationships  include  the  effects  of 
freestream  turbulence  and  streamwise  pressure  gradient.  Using  the 
relationships  of  Abu-Ghannam  and  Shaw  [1980],  for  example,  we  predicted  the 
onset  of  transition  for  the  pressure  surface  data  using  the  measured 


pressure  distribution  and  a  freestream  turbulence  intensity  (in  the  blade 
pack)  of  1.5%.  This  turbulence  intensity  value  was  determined  from  hot-wire 
measurements  at  the  edge  of  the  boundary  layer,  close  to  the  blade  leading 
edge.  Onset  of  transition  was  predicted  to  be  at  a  momentum  thickness 
Reynolds  numbers  (Reg)  of  342,  and  comparison  with  the  Reg  found  from  the 
profiles  put  this  onset  at  47.8%  chord.  By  onset  here,  we  mean  the  first 
location  at  which  the  intermittency ,  as  measured  with  a  flush-mounted  film 
probe  for  example,  would  be  greater  than  zero.  Because  of  the  strong 
favorable  gradient  near  the  trailing  edge,  the  scheme  also  predicted  that  a 
fully-turbulent  boundary  layer  would  not  develop  on  the  pressure  surface. 

There  do  not  seem  to  be  any  empirical  prediction  schemes  which  include  the 
effect  of  surface  curvature.  While  convex  curvature  apparently  has  no  effect 
on  transition,  the  concave  curvature  of  the  pressure  surface  can  promote  the 
generation  of  Gortler  vortices  (see  Gortler  [1940]),  which  can  cause 
transition  to  occur  earlier. 

Sublimation  flow  visualization  studies  helped  determine  the  transition 
point  on  the  pressure  surface.  The  average  of  five  flow  visualization  tests 
placed  the  transition  point  at  64.2 %  ±  3.9%  chord  with  95%  confidence.  Figure 
8  shows  a  plot  of  mean  velocity,  normalized  by  the  edge  velocity,  for  a  fixed 
distance  (y  =  0.508  mm)  above  the  plate  and  a  varying  chord  location.  At  this 
distance  above  the  surface,  the  measurement  volume  is  above  the  boundary  layer 
for  the  first  two  chord  locations.  The  decrease  in  mean  velocity  with  chord 
location  over  the  first  half  of  the  blade  reflects  the  growth  of  the  boundary 
layer  relative  to  the  fixed  distance.  The  rapid  rise  in  mean  velocity  near 
60%  chord  indicates  the  onset  of  transition  (see  for  example,  Klebanoff, 
Tidstrom,  and  Sargent  [1962]).  Agreement  with  the  flow  visualization  studies 
appears  to  be  quite  good.  However,  the  simple  empirically  based  calculations 


of  Abu-Ghannam  and  Shaw  [1980]  pick  too  small  a  chord  value  for  the  onset  of 
transition.  The  fact  that  the  mean  velocity  does  not  reach  a  constant  (or 
decreasing)  value  with  increasing  chord  location  indicates  that  the  transition 
process  is  not  complete. 

Integral  parameters  can  also  characterize  transitional  boundary  layers. 
Plots  of  displacement  thickness  (6*),  first  shape  factor  (Hj2)>  and  Reg  are 
shown  in  Figure  9.  Also  shown  are  values  for  the  skin  friction  coefficient 
(Cf).  Most  of  the  integral  parameters  shown  in  Figure  9  were  obtained  from 
the  smoothed  cubic  spline  approximation.  Because  of  the  lack  of  near  wall 
measurements  for  some  of  the  extremely  thin  layers,  we  felt  that  the  values  of 
momentum  thickness  obtained  from  the  spline  fit  were  not  accurate,  so  that 
some  values  of  Hj2  and  Reg  (as  shown  in  Figure  9)  were  calculated  from  the 
appropriate  Falkner-Skan  appoximations.  Note  the  large  decrease  in  6*  as  the 
flow  encounters  the  favorable  pressure  gradient  near  the  trailing  edge.  Hf2 
shows  laminar-like  values  until  just  before  the  68.0%  chord  location,  at  which 
point  the  values  drop  into  the  turbulent  regime.  The  values  of  Cf  are 
determined  from  the  Falkner-Skan  approximations.  At  the  leading  edge,  the 
skin  friction  goes  to  infinity.  In  the  transition  regime,  the  Cf  values, 
although  known  to  increase,  cannot  be  easily  estimated.  To  indicate  how  large 
Cf  might  become  near  the  trailing  edge,  a  value  based  on  the  Ludweig-Tillman 
empirical  expression  (see  Ludweig  and  Tillman  [1949])  is  given  for  the  97.9 % 
chord  location.  Use  of  the  Ludweig-Tillman  expression  here  is  not  strictly 
valid,  as  the  boundary  layer  profile  is  probably  not  fully  turbulent  at  97.9% 
chord. 

An  attempt  was  made  to  fit  the  boundary  layer  profile  at  97.9%  chord  to 
the  wall-wake  equation,  but  no  logarithmic  region  was  obtained.  Purtell, 
Klebanoff,  and  Buckley  [1981]  concluded  that  the  logarithmic  region  seems  to 
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be  an  inherent  characteristic  of  the  turbulent  boundary  layer.  That  is,  for 
fully-developed  turbulent  boundary  layers,  they  found  the  extent  of  the 

\  logarithmic  region  to  be  roughly  a  constant  fraction  of  the  boundary  layer 

i 

|  thickness  as  Reg  was  decreased.  Murlis,  Tsai,  and  Bradshaw  [1982]  found 

i 

|  strong  evidence  for  the  validity  of  the  logarithmic  law  of  the  wall,  at  zero 

pressure  gradient,  for  Reg  values  as  low  as  700,  while  Smits,  Matheson,  and 
Joubert  [1983]  found  a  logarithmic  region,  for  favorable  pressure  gradients, 
at  Reg  as  low  as  261.  Since  no  logrithmic  region  was  found  for  the  boundary 
layer  at  97.9%  chord,  it  must  be  concluded  that  either  the  boundary  layer  was 
not  fully  turbulent  at  a  Reg  of  388,  or  (less  likely)  that  the  logarithmic 
region  was  so  small  that  it  could  not  be  detected. 

Turbulence  intensity  data  for  the  pressure  surface  boundary  layer  are 
shown  in  Figure  10.  As  the  Falkner-Skan  approximation  has  been  shown  to  be 
reasonable  for  the  profiles  to  57.2%  chord,  the  large  turbulence  intensities 
near  the  wall  are  disturbing.  A  typical  profile  of  the  skewness  versus  y/6*, 
which  is  shown  for  the  5.9%  chord  location  in  Figure  11,  adds  considerably  to 
the  problem  in  that  this  profile  might  reasonably  resemble  the  shape  of  a 
skewness  profile  one  might  expect  to  find  from  measurements  of  a  turbulent 
boundary  layer. 

This  problem  was  examined  in  some  detail.  Using  both  calibrated  hot-wire 
probes  in  the  boundary  layer  and  uncalibrated  hot  films  flush-mounted  on  the 
surface,  we  determined  that  the  profiles  near  the  leading  edge  were  indeed 
laminar.  Typically,  at  the  edge  of  these  leading  edge  boundary  layers,  the 
intensity  was  found  to  be  near  1.5%.  The  difference  between  this  value  and 
the  0.18%  found  in  the  approach  flow  is  probably  due  to  the  interaction  of  the 
flow  with  the  blade  pack.  Having  shown  the  boundary  layers  to  be  laminar,  we 
next  suspected  that  the  intensity  measurements  might  be  contaminated  by 


mean-velocity-gradient  broadening.  This  problem  has  been  considered 
previously  by  Edwards,  Angus,  French,  and  Dunning  [1971],  Goldstein  and  Adrian 
[1971],  and  Kried  [1974].  For  simplicity  in  the  current  study,  the  laser 
intensity  was  taken  to  be  constant  for  the  entire  measurement  volume.  In  the 
present  case  then,  error  estimates  could  be  easily  made  by  assuming  the 
Falkner-Skan  approximate  profiles  or  by  using  the  smoothed  cubic  spline  fit. 
Similiar  results  are  obtained  for  either  estimate.  In  Figure  12a,  an  estimate 
of  the  turbulence  intensity  caused  by  mean-velocity-gradient  broadening  is 
shown  against  percent  chord  for  a  measurement  volume  roughly  as  large  as  the 
volume  estimated  in  Section  IV.  As  it  seemed  plausible  that  the  actual 
measurement  volume  might  be  larger  than  the  volume  estimated  theoretically,  we 
repeated  the  calculations  for  a  measurement  volume  roughly  twice  that  of  the 
estimated  volume.  These  results  are  shown  in  Figure  12b.  Figures  12a  and  12b 
show  estimates  for  a  constant  y/<5*  of  roughly  1.70;  the  measurement  points  are 
also  given.  It  is  clear  from  a  comparison  of  the  turbulence  intensity 
calculated  from  the  velocity-gradient  broadening  against  the  measured  data, 
that  gradient  broadening  alone  cannot  account  for  the  entire  intensity.  In 
addition,  the  skewness  when  calculated  from  an  assumption  of  velocity-gradient 
broadening  is  much  smaller  than  that  observed  experimentally. 

As  a  second  approach  to  the  problem,  we  assumed  that  in  addition  to  the 
gradient  broadening  problem,  a  small  vibration  may  have  contaminated  the 
velocity  signal.  Calculations  are  again  straightf oward  using  either  the 
spline  fit  or  the  Falkner-Skan  approximations.  Results  for  vibration 
amplitudes  of  25.4  pm  and  50.8  pm  are  again  given  in  Figure  12.  With  the 
exception  of  the  points  at  5.9%  chord,  which  appear  to  have  been  biased  by  an 
inaccurate  calculation  of  6*  (see  Figure  9),  the  calculations  agree  reasonably 
well  with  the  measurements  for  a  measurement  volume  of  66  pm  and  a  vibration 


amplitude  of  50.8  pm.  Comparison  between  the  measured  and  calculated 
skewnesses  are  also  much  closer.  Some  simple  measurements  with  an 
accelerometer  indicated  that  vibration  amplitudes  of  this  magnitude  were  not 
unreasonable  to  expect,  so  that  a  combination  of  velocity-gradient  broadening 
with  a  vibration  of  the  measurement  volume  seems  a  likely  cause  for  the 
elevated  turbulence  levels.  As  shown  in  Figures  12a  and  12b  (and  as  can  be 
shown  for  the  suction  surface  turbulent  profiles),  the  effect  becomes  quite 
small  as  the  boundary  layer  grows.  In  the  present  situation,  the  bias  can 
probably  be  considered  negligibly  small  for  chord  positions  larger  than  about 
25%. 

Turbulence  intensity  profiles  are  shown  for  the  transitioning  boundary 
layers  on  the  pressure  surface  in  Figure  10.  The  data  show  classical  shape 
(see  Klebanof,  Tidstrom,  and  Sargent  [1962])  and  agree  reasonably  well  with 
the  measurements  of  Wang,  Simon,  and  Buddhavarapu  [1985]. 

VI.  SUCTION  SURFACE  BOUNDARY  LAYERS 

Suction  surface  boundary  layers  were  taken  at  11  chordwise  locations  on 
the  center  blade.  Figure  5  shows  these  chordwise  locations  as  well  as  the 
static-pressure  distribution.  A  very  large  adverse  pressure  gradient  exists 
near  the  leading  edge.  This  gradient  gradually  becomes  less  severe  with 
downstream  distance  and  vanishes  entirely  near  80%  chord.  No  pressure 
gradient  is  evident  the  last  20%  of  the  chord  which  indicates  a  possible 
separation  region — that  is,  a  region  which  cannot  sustain  a  streamwise 
pressure  gradient. 

The  measured  suction  surface  boundary  layers  are  presented  in  Figure  13. 


As  was  the  case  for  the  pressure  surface  boundary  layers,  the  inviscid  regions 
show  the  effects  of  the  normal  pressure  gradient.  For  the  suction  surface, 
however,  3p/3y  >  0  as  opposed  to  the  pressure  surface  where  3p/3y  <  0.  Once 


again,  the  large  curvature  in  the  streamlines  near  the  leading  edge  results  in 
a  highly  curved  inviscid  region.  The  95%  confidence  bands  are  quite  small  for 
all  the  boundary  layers  except  in  two  regions.  First,  the  thin  boundary  layer 
at  2.6%  chord  has  a  large  velocity  gradient  near  the  surface;  because  of  this 
large  gradient,  the  sensitivity  to  probe  placement  is  heightened,  and  the 
measurements  are  less  repeatable.  Second,  as  suspected,  the  boundary  layer  at 
94.9%  chord  was  separated,  and  the  unsteadiness  in  the  separation  process 
resulted  in  larger  error  bands. 

All  of  the  measured  suction  surface  boundary  layers  are  turbulent.  This 
implies  that  transition  took  place  before  the  measurement  station  at  2.6% 
chord,  which  is  not  surprising  considering  the  very  large  adverse  pressure 
gradient  near  the  leading  edge.  The  separation  of  a  laminar  boundary  layer 
under  an  adverse  pressure  gradient  results  in  a  free  shear  layer,  which  is 
unstable.  The  transition  to  turbulence  takes  place  very  rapidly.  Once 
turbulent  entrainment  increases,  the  shear  layer  is  enlarged  which  results  in 
a  pressure  recovery  and  a  rapid  reattachment.  Thus,  the  separation  "bubble" 
can  be  quite  short  and  close  to  the  leading  edge. 

Although  the  transition  takes  place  very  close  to  the  leading  edge,  the 
recovery  process  extends  some  distance  downstream.  This  process  can  be  seen 
from  the  mean-velocity  profiles  plotted  in  dimensionless  outer  variables. 

These  plots  are  shown  in  Figure  14  where  the  normal  pressure  gradients  have 
been  taken  into  account  as  described  previously  (see  Zierke  and  Deutsch 
[1985]).  The  recovery  process  can  be  seen  to  extend  through  the  2.6%  and  7.6% 
chord  locations  by  observing  the  shape  of  the  profiles.  As  we  will  show 
later,  the  shape  of  the  velocity  profiles  results  in  higher  values  of  Hj2 
during  the  recovery  process.  The  local  turbulence  intensity  profile  at  2.6% 
chord  also  indicates  recovery.  This  turbulent  boundary  layer  was  the  only  one 


measured  in  which  there  was  a  maximum  value  of  local  turbulence  intensity  away 
from  the  surface. 

The  mean  velocity  data  were  fit  to  the  wall-wake  equation  of  Coles 
[1956], 


through  a  least  squares  technique  described  by  Zierke  and  Deutsch  [1985].  For 
a  given  boundary  layer  thickness,  the  technique  simultaneously  calculates  the 
values  of  uT  and  II  which  yield  the  best  fit  to  the  data.  Figure  15  shows  the 
velocity  profiles  in  inner  variables.  The  logarithmic  region  reaches  a 
maximum  and  the  wake  region  reaches  a  minimum  at  12.7%  chord.  This  seems  to 
be  a  second  indication  of  complete  recovery  from  the  leading  edge  separation 
“bubble."  Moving  further  downstream.  Coles'  wake  parameter,  II,  (which 
controls  the  size  of  the  wake  region)  and  Reg  increase  resulting  in  a 
reduction  in  the  extent  of  the  logarithmic  region.  As  separation  is  reached, 
the  logarithmic  region  disappears  and  the  wall-wake  equation  cannot  be  fit  to 
the  data.  This  conclusion  was  reached  earlier  by  Simpson,  Strickland,  and 
Barr  [1977],  who  found  the  law  of  the  wall  valid  until  intermittent  separation 
(flows  containing  instantaneous  flow  reversals)  was  reached. 

Although  the  influence  of  surface  curvature  cannot  be  extracted  from  the 
data,  one  must  suspect  that  this  influence  is  indeed  present.  The  convex 
curvature  on  the  suction  surface  (0.01  <  |6/RC|  <  0.2)  and  the  concave 
curvature  on  the  pressure  surface  (0.002  <  |  61^.  (  <  0.02)  have  opposite  effects 
on  turbulent  boundary  layers.  Ramaprian  and  Shivaprasad  (1977]  show  that 
convex  curvature  reduces  the  logarithmic  region  and  increases  the  relative 
strength  of  the  wake  component.  Except  for  the  initial  region  of  curvature, 
convex  curvature  Increases  the  rate  of  growth  of  Reg  and  decreases  Cf. 
Shivaprasad  and  Ramaprian  [1978]  claim  that  the  effects  of  convex  curvature  on 
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the  behavior  of  the  turbulent  boundary  layer  are  even  stronger  than  the 
effects  of  concave  curvature  at  the  same  value  of  [6/Rc|.  Their  measurements 
showed  that  convex  curvature  reduces  turbulence  intensity  and  Reynolds  shear 
stress.  Measurements  by  So  and  Mellor  [1973],  Gillis  and  Johnston  [1983],  and 
Gibson,  Verriopoulos,  and  Vlachos  [1984]  agreed.  These  results  indicate  that 
for  very  strong  convex  curvature  effects,  regions  can  be  found  where 
turbulence  cannot  exist.  Bradshaw  [1969]  showed  that  the  behavior  of  the 
turbulent  boundary  layer  is  very  sensitive  to  streamline  curvature  as  mild  as 
|6/RC|  =  0.003.  He  used  an  analogy  between  the  effects  of  streamline 
curvature  and  buoyancy  to  estimate  quantitatively  the  effect  of  curvature  on 
mixing  length  distribution  in  the  boundary  layer.  So  [1975]  verified  this 
buoyancy  analogy  mathematically.  Shivaprasad  and  Raraaprian  [1978]  made 
measurements  which  support  the  buoyancy  analogy  of  Bradshaw  [1969]  for  mild 
convex  curvature.  For  concave  curvature,  they  found  the  bouyancy  analogy 
useful  only  for  values  of  | 6/Rc |  near  0.01. 

Figure  16  includes  plots  of  6*,  Hj2,  Re©,  and  Cf.  These  parameters  were 
calculated  from  a  smoothed  cubic  spline  fit  of  the  data,  except  for  Cf,  which 
was  calculated  from  the  least  squares  fit  of  the  data  to  the  wall-wake 
equation.  Values  of  and  Re©  were  also  used  to  calculate  Cf  from  the 
empirical  equation  of  Ludweig  and  Tillman  [1949].  The  displacement  thickness 
increases  gradually  at  first  and  then  increases  rapidly  through  separation. 

The  plot  of  Hj2  indicates  a  turbulent  boundary  layer  beginning  near  the 
leading  edge.  Recovery  from  the  leading  t dge  separation  "bubble”  results  in 
an  initial  decrease  of  H j 2 •  Separation  of  turbulent  boundary  layers  is 
usually  approximated  using  values  of  H j 2  near  2.2  which  corresponds  here  to  a 
suction  surface  location  near  60%  chord.  Sandborn  and  Kline  [1961]  proposed  a 
relation  for  intermittent  separation. 
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which  yields  66.9%  chord  (corresponding  to  =  2.70)  as  the  location  of 
intermittent  separation  for  the  data  presented  here.  Values  of  Cf  appear  to 
be  near  zero  at  the  leading  edge  which  corresponds  to  the  vanishing  skin 
friction  at  the  beginning  of  the  leading  edge  separation  "bubble."  Cf  reaches 
a  maximum  after  the  boundary  layer  has  totally  recovered  from  the  leading  edge 
separation  and  then  decreases  as  the  trailing  edge  separation  of  the  turbulent 
boundary  layer  is  reached.  Cf  vanishes  near  80%  chord. 

Defining  separation  as  the  entire  process  of  the  breakdown  of  boundary 
layer  flow,  Simpson,  Chew,  and  Shivaprasad  [1981]  quantified  the  various 
stages  of  separation  with  the  instantaneous  backflow  near  the  wall.  Incipient 
detachment  (ID)  occurs  with  1%  instantaneous  backflow;  intermittent  transitory 
detachment  (ITD)  occurs  with  20%  instantaneous  backflow;  transitory  detachment 
(TD)  occurs  with  50%  instantaneous  backflow;  and  detachment  occurs  when  the 
wall  shear  stress  becomes  zero.  The  percent  backflow  is  easily  calculated  as 
the  portion  of  the  measured  velocity  distribution  that  includes  negative 
velocities.  Figure  17  shows  the  instantaneous  backflow  measurements  at  the 
63.2%,  74.0%,  84.2%,  and  94.9%  chord  locations.  Figure  18  shows  the  maximum 
percent  backflow  as  a  function  of  percent  chord.  The  sublimation  flow 
visualization  tests  showed  separation  to  occur  at  65.6%  ±  3.5%  chord,  and  a 
comparison  with  the  maximum  instantaneous  backflow  data  of  Figure  18  shows 
that  flow  visualization  yields  a  value  for  separation  which  is  only  slightly 
downstream  of  incipient  detachment.  Locating  turbulent  separation  by 
observing  when  Hj2  nears  2.2  also  seems  to  indicate  incipient  detachment.  We 
might  note  that  although  Simpson,  Chew,  and  Shivprasad  [1981]  state  that 
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detachment  and  transitory  detachment  need  not  be  at  the  same  location,  our 
skin  friction  calculations  show  that  the  chordwise  locations  of  detachment  and 
transitory  detachment  are  quite  close  to  one  another. 

Restrictions  in  applying  the  wall-wake  equation  in  the  vicinity  of 
separation  result  from  the  velocity  scale,  uT,  approaching  zero.  A  vanishing 
uT  leads  to  a  vanishing  logarithmic  region,  which  would  not  cause  concern  if 
the  separated  flow  profiles  would  follow  the  law  of  the  wake.  Unfortunately, 
experimental  data  prove  otherwise  (see  Simpson,  Chew,  and  Shivaprasad  [1981] 
for  example).  Perry  and  Schofield  [1973]  developed  a  similarity  defect  law 
based  on  the  maximum  shear  stress  rather  than  the  wall  shear  stress.  The 
defect  law  was  originally  developed  for  attached  boundary  layers  under 
moderate  to  strong  adverse  pressure  gradients  where  xmax/rw  >  1.5.  The  defect 
law  is 


where 

Us  =  8.0 

and 

B  -  2.86  6*  . 

Um  is  a  velocity  scale  based  on  the  maximum  shear  stress  (Um  =  ^imax/P )  and  L 
is  the  distance  from  the  wall  to  Tmax.  The  velocity  scale  Us  is  found  using 
the  methodology  that  Clauser  [1954]  used  to  determine  uT.  A  half-power 
equation  is  used  near  the  wall, 
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Perry  and  Schofield  [1973]  suggested  using  the  defect  law  for  the  outer  90%  of 
the  boundary  layer  with  the  half-power  equation  forming  the  innermost  portion 
of  that  defect  law.  They  recommended  the  law  of  the  wall  as  an  inner  wall 
matching  condition.  Schofield  [1983]  extended  the  model  to  detached  flows  by 
suggesting  that  the  similarity  would  hold  provided  that  the  origin  of  the 
normal  coordinate  has  been  moved  from  the  wall  to  the  location  at  which 
u  =  0. 

The  suction  surface  boundary  layers  are  plotted  using  the  similarity 
relation  of  Perry  and  Schofield  [1973]  in  Figure  19.  In  the  outer  90%  of  the 
boundary  layer,  the  similarity  relation  collapses  the  data  quite  well  for 
chord  locations  upstream  of  the  separated  region.  As  the  amount  of 
instantaneous  backflow  increases,  however,  the  data  deviates  more  and  more 
from  the  similarity  relation.  This  deviation  is  in  seeming  contrast  to  the 
conclusions  made  by  Schofield  [1983].  A  close  examination  of  his  defect 
plots,  however,  shows  similar  trends  in  his  analyzed  data  and  the  data  shown 
in  Figure  19. 

No  outer  region  similarity  seems  to  exist  downstream  of  detachment.  Many 
researchers  have  attempted  a  law  of  the  wake  correlation  without  success  (see 
Simpson,  Chew,  and  Shivaprasad  [1981]  for  example).  For  a  turbulent  boundary 
layer  subjected  to  a  streamwise  pressure  gradient,  Mellor  and  Gibson  [1969] 
suggested  replacing  the  shear  velocity  with  a  pressure  velocity, 

( 6*/p)(dp/dx),  as  the  velocity  scale.  However,  Schofield  [1981]  has  shown 
this  scaling  to  be  inadequate.  Using  the  velocity  scale  corresponding  to  the 
maximum  shear  stress  as  suggested  by  Schofield  [1983]  has  been  shown  to  give 
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1  poor  similarity  in  Figure  18.  Mehta  and  Goradia  [1984]  had  some  success  by 

t 

j  assuming  the  the  outer  region  velocity  profiles  behave  like  a  two-dimensional 

mixing  layer.  Their  similarity  variables  were  not  found  to  give  outer  region 
similarity  with  the  data  measured  here. 

|  Similarity  in  the  backflow  region  seems  to  show  more  promise.  Simpson, 

Chew,  and  Shavaprasad  [1981]  found  good  backflow  similarity  by  normalizing  the 
velocity  by  the  maximum  backflow  velocity,  and  the  distance  from  the  wall  by 
the  distance  to  the  maximum  backflow  velocity.  Schofield  [1983]  found  that 
this  backflow  similarity  could  be  improved  by  using  the  total  backflow 
thickness  as  the  length  scale.  Figure  20  shows  this  backflow  similarity  for 
the  data  at  94.9%  chord.  Despite  the  scatter,  the  backflow  data  seems  to 
collapse  quite  well  with  the  data  measured  by  Simpson,  Strickland,  and  Barr 
[1977]  and  Simpson,  Chew,  and  Shivaprasad  [1981].  The  only  exceptions  are  the 
two  data  points  closest  to  the  wall — data  for  which  the  95%  confidence  bands 
are  larger  than  the  magnitude  of  the  mean  velocity. 

The  turbulence  intensity  on  the  suction  surface  is  shown  for  all  eleven 
chord  positions  in  Figure  21.  Recovery  from  the  leading  edge  separation  is 
apparent  from  the  peaks  in  turbulence  intensity  that  occur  away  from  the 
surface  for  the  2.6%  and  7.6%  chord  locations.  These  peaks  also  occur  in  the 
separation  region  from  63.2%  chord  through  94.9%  chord.  This  shape  reflects 
the  movement  In  the  location  of  the  maximum  mean-shear  rate  outward  from  the 
near-wall  region. 

vn.  WAKES 

Near-wake  measurements  were  made  at  105.4%  chord  and  109.6%  chord  using 
the  LDV  technique.  Five-hole  probes  were  used  to  measure  the  far  wake  at 
152.6%  chord.  Figure  22  shows  the  data  points  and  their  95%  confidence  bands 
for  all  three  wake  profiles.  The  two  near-wake  profiles,  which  are  quite 
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j  similar,  are  very  asymmetric  because  of  the  large  difference  in  trailing  edge 

|  boundary  layer  thicknesses  on  the  two  surfaces  of  the  blade.  The  separation 

of  the  suction  surface  boundary  layer  leads  to  negative  mean  velocities  at  the 
center  of  the  near  wake.  Other  researchers  have  also  measured  negative  mean 
velocities  in  near  wakes.  Wadcock  [1980],  using  a  flying  hot  wire,  measured 
negative  mean  velocities  in  the  near  wake  of  an  airfoil.  Braden,  Whipkey, 
Jones,  and  Lilley  [1983]  used  a  LDV  to  measure  negative  mean  velocites  in  the 
near  wake  of  an  airfoil  with  confluent  boundary  layers.* 

Wakes  become  similar  only  at  distances  far  downstream  of  their  source.  A 
Gaussian  distribution  can  be  used  to  correlate  these  far-wake  data. 
Lakshminarayana  and  Davino  [1979]  suggested  the  correlation 

ue  ~  u  -0.693  n 
Ue  ~  UCL 

where  n  is  the  normalized  distance  across  the  wake.  The  suction  and  pressure 
sides  of  the  wake  use  different  length  scales,  Ls  and  Lp.  Ls  and  Lp  are  the 
distances  on  the  suction  and  pressure  sides  of  the  wake  centerline  from  the 
point  of  minimum  velocity  to  a  point  where  the  velocity  defect  is 
(Ue  -  ucl)/2.  The  far-wake  data  of  Lakshminarayana  and  Davino  [1980]  showed 
similarity  away  from  the  wake  edge  for  both  inlet  guide  vane  wakes  and  stator 
blade  wakes.  This  similarity  was  corroborated  by  Hobbs,  Wagner,  Dannenhof fer, 
and  Dring  [1982]  for  distances  greater  than  30%  axial  chord  downstream  of  the 
trailing  edge  of  their  compressor  cascade  blades.  Figure  23  shows  that  our 
wake  data  has  Gaussian  similarity  in  the  far  wake.  However,  the  two  near-wake 
profiles  exhibit  no  similarity. 

^Confluent  boundary  layers  develop  on  an  airfoil  with  leading  edge  slats  or 
trailing  edge  flaps  which  causes  the  boundary  layers  from  the  various 
surfaces  to  interact. 
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Integral  parameters  can  be  calculated  from  the  wake  velocity  profiles. 
Using  the  location  of  the  wake  centerline  and  displacement  thicknesses  on 
either  side  of  the  centerline,  the  displacement  surface  can  be  drawn  as  seen 
in  Figure  24.  The  curvature  of  the  displacement  surface  in  the  near-wake 
region  can  be  explained  by  the  pressure  difference  of  the  two  blade  surfaces. 
The  location  of  the  displacement  surface  in  the  far-wake  region  can  be  partly 
explained  by  the  large  amount  of  separation  on  the  suction  surface  of  the 
neighboring  blades.  However,  the  major  explanation  for  this  curvature  of  the 
displacement  surface  in  the  far-wake  region  is  the  confinment  of  the  far  wakes 
between  the  two  tailboards  (see  Figure  2).  The  effects  of  the  tailboards  must 
be  taken  into  account  when  considering  this  displacement  surface  for  wake 


modeling. 

Turbulence  intensity  profiles  are  presented  for  the  two  near  wakes  in 
Figure  25.  As  with  the  turbulence  intensity  profiles  in  the  separating 
boundary  layers,  the  turbulence  intensity  peaks  are  displaced  outward, 
essentially  tracking  the  regions  of  large  mean-velocity  gradients.  These 
data,  although  more  detailed,  are  quite  similar  to  the  data  of  Hah  and 
Lakshminarayana  [1982]  for  the  near  wake  of  an  isolated  airfoil. 


VIII.  CONCLUSIONS 

Viscous  calculations  for  turbomachinery  applications  have  been 
handicapped  by  a  lack  of  sufficiently  detailed  and  precise  data  against  which 
these  calculation  schemes  can  be  compared.  In  order  to  help  overcome  this 
problem,  we  have  presented  measurements  of  the  boundary  layers  and  wakes  about 
a  double  circular  arc,  compressor  blade  in  cascade.  A  two-dimensional, 
periodic  cascade  flow  has  been  developed  without  the  use  of  continuous  side 


wall  suction.  This  facility  has  allowed  these  measurements  to  be  made  with  a 
non-intrusive  LDV  system. 
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Despite  the  facts  that  the  measured  flow  field  was  very  complex  and  that 
only  a  one-component  LDV  system  was  used,  the  measurements  presented  here  have 
led  to  an  initial  physical  understanding  of  the  cascade  flow  field.  Regions 
of  laminar  flow,  transition,  recovery  from  a  leading  edge  separation  "bubble," 
non-equilibrium  turbulent  flow,  separated  flow,  and  near-wake  flow  have  all 
been  investigated.  Prediction  of  a  flow  field  that  includes  all  of  these  flow 
regions  will  be  a  severe  test  for  any  viscous  computational  technique. 

Despite  the  varying  streamwise  pressure  gradient,  the  laminar  velocity 
profiles  near  the  leading  edge  on  the  pressure  surface  show  reasonable 
agreement  with  Falkner-Skan  velocity  profiles  (computed  at  the  same  streamwise 
pressure  gradient).  Transition  was  identified  through  a  departure  of  the 
measured  profile  shape  from  the  Falkner-Skan  profile  shape,  or  through  the 
boundary  layer  shape  factors.  Sublimation  flow  visualization  tests  agree  well 
with  the  position  of  the  transition  region,  but  empirical  relationships 
predict  onset  of  transition  somewhat  early.  Transition  on  the  pressure 
surface  was  incomplete. 

A  problem  was  encountered  for  the  LDV  measurements  in  the  extremely  thin 
laminar  boundary  layers  on  the  pressure  surface.  Large  turbulence  intensities 
were  indicated.  This  problem  was  traced  to  a  combination  of  mean-velocity- 
gradient  broadening  and  measurement  volume  vibration.  Turbulence  intensity 
profiles  in  the  transitional  and  turbulent  boundary  layers  were  not  affected. 

The  leading  edge  separation  "bubble"  on  the  suction  surface  was  too  small 
to  be  measured.  Boundary  layers  measured  downstream  of  this  "bubble"  are 
fully  turbulent  and  the  recovery  process  after  reattachment  extends  downstream 
a  distance  of  nearly  10%  chord.  The  recovery  process  can  be  identified  by  the 
velocity  profile  shapes  (including  shape  factors). 


The  non-equlibrium  turbulent  boundary  layers  downstream  of  the  separation 
bubble"  follow  the  wall-wake  equation  of  Coles  [1956]  until  detachment  of  the 


boundary  layer  is  reached.  These  boundary  layers  also  show  good  similarity 
using  the  defect  law  of  Perry  and  Schofield  [1973]  up  to  the  location  where 
Instantaneous  backflow  is  present.  Similarity  becomes  worse  as  the  amount  of 
instantaneous  backflow  Is  Increased.  No  outer  region  similarity  seems  to 
exist  downstream  of  detachment.  In  the  backflow  region,  however,  the  data 
measured  here  seem  to  follow  the  backflow  similarity  shown  by  other 
researchers.  Backflow  similarity  is  found  by  using  the  maximum  backflow 
velocity  as  the  velocity  scale  and  the  total  backflow  thickness  as  the  length 
scale. 

The  near-wake  velocity  profiles  are  asymmetric  and  include  negative  mean 
velocities  at  the  wake  center.  These  profiles  do  not  show  the  Gaussian 
similarity  shown  in  the  far-wake  profile  measured  with  a  five-hole  probe. 
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X.  APPENDIX 


Equations  for  the  Double  Circular  Arc  Blades 


The  equations  for  the  pressure  surface,  suction  surface,  and  camberline 
of  the  double  circular  arc  blades  used  in  the  current  study  can  be  written  as 

2  2  2 
XP  +  [yp  +  219.7]  =  246.8 


2  r  -i 2  2 

xs  +  [ys  +  149.5]  =  189.1  , 


2  r  -i 2  2 

xc  +  [yc  +  179.4]  =  212.8 


All  of  the  dimensions  are  in  millimeters.  The  origin  of  the  coordinate  system 
used  here  is  located  on  the  chord  line  of  midchord.  The  x-coordinate  is 
parallel  to  the  chord,  while  the  y-coordinate  is  normal  to  the  chord. 
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Nomenclature 


A 

B 

c 

C 

Cf 

CP 

D 

H12 

i 

ID 

ITD 

L 

LP>  l^s 


LDV 

M 

n 

N 

P 

?T 

r 

Rc 

Rec 

Re0 

s 

SFV 


aspect  ratio 

Integral  layer  thickness  In  the  Perry  and  Schofield  [1973]  theory 

blade  chord  length 

law  of  the  wall  constant  (=  5.0) 

skin  friction  coefficient  *  Tw/(p  Ue^/2) 

static  pressure  coefficient  *  (p  -  pi)/(p  V\^/2) 

diffusion  factor 

first  shape  factor  =  6*/0 

incidence  angle  =  gj  - 

incipient  detachment 

intermittent  transitory  detachment 

distance  from  the  surface  to  the  location  of 

pressure  and  suction  surface  length  scales  from  the  point  of 
maximum  velocity  to  a  point  where  the  velocity  detect  is 
(Ue  -  uCL)/2 

laser  Doppler  velocimeter 

Mach  number 

data  point  index 

number  of  data  points 

static  pressure 

total  or  stagnation  pressure 

radius 

radius  of  curvature 
olade  chord  Reynolds  number  =  cVj/v 
momentum  thickness  Reynolds  number  =  6Ue/v 
blade  spacing 

location  of  separation  from  flow  visualization  tests 


Nomenclature  (Cont  'd  ) 


transitory  detachment 
streamwise  velocity 

root  mean  square  value  of  the  turbulent  velocity  fluctuation 
dimensionless  velocity  in  the  inner  boundary  layer  ■  u/uT 
shear  or  friction  velocity  ■  /tw/ p 
maximum  backflow  velocity 

velocity  at  the  boundary  layer  or  wake  edge 
velocity  scale  based  on  the  maximum  shear  stress  =  /xraax/p 
velocity  scale  for  the  Perry  and  Schofield  [1973]  defect  law 
velocity 

Coles'  universal  wake  function  =  2  sin^(  )  =  1  -  cos(  ) 

streamwise  coordinate;  blade  coordinate  (see  Appendix) 

coordinate  normal  to  the  blade  surface  or  across  the  wake;  or 
blade  coordinate  (see  Appendix) 

dimensionless  coordinate  normal  to  the  blade  surface  in  the 
inner  boundary  layer  =  yuT/v 

flow  angle  measured  from  the  axial  direction 

stagger  angle 

boundary  layer  thickness  (where  u  =  0.99  Ue) 


displacement  thickness 


/  (1  -?-)  dy 

o  Ue 


depth  of  backflow 
deviation  angle  =  &2  ~  K2 

fluid  turning  or  deflection  angle  =  fcS ^  —  02 


normalized  distance  across  the  wake 


Nomenclature  (Cont  'd ) 


> 

CO 

momentum  thickness  =  /  yy—  ( 1  -  yy— )  dy 

o  e  e 

c 

von  Karman 's  mixing  length  parameter  (*  0.41);  blade  metal  angle 

\> 

kinematic  viscosity  (0.150  cm2/sec  for 

air) 

1 

Coles'  wake  parameter 

p 

fluid  density  (1.205  kg/m3  for  air) 

a 

blade  solidity  =  c/s 

Tmax 

maximum  shear  stress 

Tw 

wall  or  surface  shear  stress 

<t> 

camber  angle  *  <1  -  <2 

u 

total-pressure  loss  coefficient  =  (P^ 

-  P t  )/ (p 

Vi2/2) 

Subscrij 

>ts 

c 

camber  line 

CL 

at  the  wake  centerline 

inv 

inviscid 

LE 

leading  edge 

m 

mean  flow 

meas 

measured 

n 

data  point  index 

P 

pressure  surface 

s 

suction  surface 

TE 

trailing  edge 

X 

axial  direction 

e 

tangential  direction 

Nomenclature  (Cont'd) 


inlet  (upstream  five-hole  probe  measurement  station) 
outlet  (downstream  five-hole  probe  measurement  station) 


Superscript 


average  over  the  blade  passage 


V’J 
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Figure  2.  Cascade  Test  Section  with  Flow  Controls 


Figure  3.  Five-Hole  Probe  Outlet  Survey 


Measured  Pressure  Surface  Boundary  Layers 
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Figure  10. 


Turbulence  Intensity  Data  for  the  Pressure  Surface 
Boundary  Layers 
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Figure  12.  Pressure  Surface  Turbulence  Intensities  at  y/6*  -  1.70,  Both 
Measured  Data  and  Data  Estimated  by  Mean-Velocity-Gradient 
Broadening  and  Measurement  Volume  Vibration  (  -  Measured 

Data;  -  Vibration  Amplitude  of  0  pm;  -  Vibration  Amplitude 
of  25.4  ym;  Vibration  Amplitude  of  50.8  ym) 
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Figure  20.  Backflow  Similarity  of  the  Reconstructed  Boundary  Layer 
at  94.9%  Chord  on  the  Suction  Surface  (The  Solid  Line 
Represents  the  Backflow  Similarity  of  Data  Measured  by 
Simpson,  Strickland,  and  Barr  [1977]  and  Simpson,  Chew 
and  Shivaprasad  [1981]) 
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Figure  21.  Turbulence  Intensity  Data  for  the  Suction  Surface  Boundary 


